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Abstract 

We show that in d > 2 dimensions the A^-particle kinetic energy 
operator with periodic boundary conditions has symmetric eigenfunc- 
tions which vanish at particle encounters, and give a full description of 
these functions. In two and three dimensions they represent common 
eigenstates of bosonic Hamiltonians with any kind of contact inter- 
actions, and illustrate a partial 'multi-dimensional Bethe Ansatz' or 
'quantum KAM theorem'. The lattice analogs of these functions exist 
for N < Lt"^/^! where L is the linear size of the box, and are common 
eigenstates of Bose-Hubbard Hamiltonians and spin-^ XXZ Heisen- 
berg models. 
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1 Introduction 



In this paper we construct eigenstates to Hamiltonians of the form 

N 

Hn = -5^Ai + V(ri,...,rjv) -/xA^ 

1=1 

V{ri, . . . , tn) = if all are different, (1) 

that is, Hamiltonians with real symmetric contact interactions. The defini- 
tion of such operators with the use of Dirac deltas is straightforward in one 
dimension and leads to integrable systems [1] . In higher dimensions the mere 
definition becomes mathematically tricky. The two- and three-dimensional 
realization of (1) needs renormalization to zero of the interaction strength 
[2], while in four dimensions and above the family (1) reduces to a unique 
element with V = 0, i.e. the kinetic energy operator ([3], Theorem X.ll). 
One can gain an impression about how this conclusion emerges by observing 
that the existence of contact interactions is closely related to the existence 
of a Green's function for the Laplacian, i.e. an L^-solution Gz{- — r') of the 
inhomogeneous equation 

(-Ar-z)G.(r-r')=5(r-r') (2) 

for any complex (non-real) z. Any Hamiltonian with a contact interaction is 
a self-adjoint extension of the N-particle Laplacian defined on C°°-functions 
whose support avoids the set {rj — rj,i j}. This symmetric operator must 
have non-vanishing deficiency indices, i.e. its adjoint must have eigenstates 
with complex eigenvalues. For two particles, these eigenstates are the Green's 
functions for complex values of z. If an L^-solution exists for equation (2), 
it has a Fourier representation 

1 r „ip{r-r') 

In one dimension this integral is convergent for any r— r' and yields a bounded 
continuous function. In two and three dimensions J |p^ — ^|~^dp is finite. 
Therefore, by Plancherel's theorem, (p^ — z)~^ is the Fourier transform of an 
L^-function which, however, diverges at r = r' as In |r — r'| and |r — r'|~^, 
respectively. In four dimensions and above / |p^ — ;i;|~^dp = oo and, thus, (2) 
has no square-integrable solution. Now by separating off the motion of the 



2 



center of mass, the two-body problem with a (5-interaction can be reduced to 
a one-body problem in a potential a5. Multiplication of a suitable function 
/ with a5 has to be interpreted as a (5, f)5, where (5, ■) denotes the action 
of the linear functional 5. If we want to define A = —A + a{6, ■)6 as a 
self-adjoint operator on L^(]R'^), we have to restrict it to functions t/j E L'^ 
such that Alp is also in and {ip, Aip) is real. So Aip cannot contain a S, 
and this is why has to be included in ip. Actually z = i suffices, and a 
straightforward computation yields that a general ip in the domain of A is of 
the form 

V'(r) = 7p{r) + a^PiO) Re Gi{r) (4) 

with 

" l + a{d,ReGi) ^ ' 

and then 

Alp = -A'ip + ai){'d)ReiGi . (6) 

Here ip is any vector from the domain of the kinetic energy —A. The problem 
appears in the denominator of a, with (5, ReGj). This is well-defined in one 
dimension and yields ReGi(O) = In dimensions two and three 5 is not 
defined on ReG^. One can proceed at least in two different ways. One may 
define 5 so that it does not feel the singularity of Gi and set (5, Re Gi) — c. 
Then different constants yield different operators, e.g c = oo gives A = —A. 
This more recent method was followed in [4]. The original procedure of [2] 
replaced (5 by a sequence Sn of regular approximants and a by When 
{^n, ReGi) — > oo. In order to obtain a nonvanishing limit for a, an 
then has to go to zero; otherwise we get back the kinetic energy operator. 

The general, A^-particle case is much more involved since it requires con- 
sidering multiple collisions and operators with infinite deficiency indices. The 
crucial question is to find the physically meaningful (lower semibounded) self- 
adjoint extensions. For this purpose generalized contact interactions have to 
be introduced [5] , [7] . For an exhaustive mathematical treatment the reader 
can consult the monographs [6] and [7] which also contain a detailed bibliog- 
raphy of the abundant mathematical and physical literature. As an example 
from the latter, we mention the use of contact interactions as 'Fermi pseu- 
dopotentials' for an approximate treatment of the hard-core Bose gas [8]. 
They also appear in connection with the Laughlin wavcfunction [9] and the 
trapped boson gas [10], [11], in which case the ground state of the projected 
(to the lowest oscillator level) Hamiltonian for fixed nonvanishing angular 
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momenta was found. Projection to the lowest oscillator level has a regular- 
izing effect on the interaction, which is however not sufficient to settle all 
problems of definition. To avoid the lower unboundedness of the projected 
Hamiltonian and to make sure that the ground state energy does not decrease 
with an increasing angular momentum, the interaction strength still has to 
decrease as 1/N with the number of particles, see equations (4) and (5) of 
Ref. [fO]. We note that apparently all existing mathematical results concern 
systems of particles in infinite space, and cannot but anticipate the extension 
of these results to confined systems. 

We are interested in the construction of eigenstates to the Hamiltonians 
(1) for the following reason. The delta-gas in one dimension is an integrable 
system, and we would like to see whether in two and three dimensions it 
preserves some trace of integrability. A partial survival of integrability could 
be viewed either as a remnant of the Bethe Ansatz when the dimensionality 
is increased or as a germ of a quantum-KAM theorem when, in the free Bose- 
gas, a contact interaction is switched on. The Bethe Ansatz solution of the 
one-dimensional delta-gas assigns N wave numbers to each eigenf unction. 
The two and three dimensional analogy would be the existence of eigenstates 
which can be expressed in terms of finitely many wave vectors. From the 
point of view of the KAM theory, surviving, maybe slightly distorted, free- 
particle eigenstates would correspond to classical KAM tori [12]. What we 
actually find are eigenstates of N bosons in a box of linear size L which 
remain unperturbed by contact interactions, are characterized by wave 
vectors and have an energy at least of order N^/L"^. 

A common feature of the operators (1) is that, on smooth functions which 
vanish at particle encounters, all act as the pure kinetic energy operator. 
Therefore, if this latter has eigenstates vanishing whenever = for some 
i 7^ j, these will be common eigenstates of all members of the family (1). For 
two particles in two dimensions it is very easy to give examples. Probably 
the simplest of them are 

2mT . . 2m7T . 2n7r , . 2n7r 

sm^-{xi-X2) sm^—{yi-y2) and cos -^(xi-X2)-cos -^(l/i-y2) ■ 

We will show that such eigenstates exist for any > 2 in any dimension 
d>2 and will find all of them for a system of bosons confined in a L x • • • x L 
torus A. That is, we will describe all A^-particle eigenstates t/j{ri, . . . ^r^) of 
the d > 2-dimensional kinetic energy operator with the following properties: 
(PI) ip is L-periodic separately in each of the dN coordinates, (P2) ip is 
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a symmetric function of ri, ... , r^v and (P3) . . , r, . . . , r, . . .) = 0. To 
satisfy (P2) and (P3) at the same time, a very large number of degenerate 
free-particle eigenstates will have to be combined. This becomes possible 
because the kinetic energy operator is unbounded and so is the degeneracy 
of the eigensubspaces. On the other hand, this explains why the energy of 
these states must increase so fast with N. As the proof of Theorem 2 will 
show, an interesting interpretation can be given to these functions: they can 
be considered as stationary states of N impenetrable bosons on the rf-torus 
which interact with each other through elastic collisions, see equations (19) 
and (25) below. If -01,-02, •■ ■ is an orthonormal basis for (Pl-3) eigenstates 
with corresponding energy eigenvalues Ei < E2 < ■ ■ ■, we can define a self- 
adjoint operator 

00 

-f^elastic = Ei\tjji) {ipi I (7) 
i=l 

which, hence, is associated to such a system. Also, every member of the 
family (1) can be decomposed as Hn = iJ^ + i/ciastic, with H'j^ depending on 
the interaction and living on the orthogonal complement of the subspace of 
(Pl-3) eigenstates. 

The problem in continuous space, that we have discussed above, will be 
presented in Section 2. In Section 3 we separately treat the lattice case. 
We shall consider only hypercubic lattices. The analog of (1) is obtained 
by interpreting A, as the lattice Laplacian. In lattices contact interactions 
mean on-site interactions which exist in any dimension and define the so- 
called Bose-Hubbard Hamiltonians. Incidentally, we shall obtain common 
exact eigenstates in any dimension d > 2 for a larger family of Hamiltonians, 

-f^lattice = ± Xl("r' " a*r){ar' " ^r) + ^ /r(nr) - /X ^ Ur 
(rr') r r 

d 

+ X] /r,r'(^r,?lr') ■ (8) 

i=l (r-r')llei 

In the first term the summation goes over nearest neighbor pairs, a* and a,. 
create, resp., annihilate a boson at site r, rir = a*ar, is the unit vector 
along the ith coordinate axis, /x is the chemical potential and fr and fr,r' 
are real functions with fr{n) = if n = 0, 1 and fr,r'{m,n) = if mn — 
0. Among others, every spin-| nearest-neighbor XXZ Heisenberg model, 
including the planar and the ferro- and antiferromagnetic isotropic models. 
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is a member of this family (with fj.{t) — +00 ii t > 2). However, there are 
two Umitations in the lattice case that we do not meet in the continuous 
one. First, the condition for finding the eigenstates is a set of transcendental 
Diophantine equations, in contrast with the algebraic ones in the continuous 
space version. We find some solutions but cannot pretend to find all. Second, 
the solutions we find are only for N < Ll'^/^l where [•] means integer part, and 
this seems to be an intrinsic limitation, due to the boundedness of i^iattice- 
In hard-core models, particle-hole transformation yields eigenstates also for 
AT > |A| - LW^I 

Finally, we note that both in the continuum and on lattices there is a 
family of exact eigenstates whose energy tends to zero as the volume in- 
creases. In two and three dimensions these states may contain N = o{L'^^^) 
particles. They correspond to gapless excitations above the grand canonical 
ground state if this latter belongs to = and, thus, has zero energy. This 
is the case e.g. if the interaction is repulsive and the chemical potential is 
zero or negative (and we take the plus sign in (8)), but may also occur with 
attractive interactions, as in the isotropic Heisenberg ferromagnet. 

2 Bosons in continuous space 
2.1 Theorems and proofs 

For spinless fermions all states share the property (P3). This observation is 
the starting point of the construction of such eigenstates for (spinless) bosons 
in more than one dimension. 

We introduce the short-hands R = (ri, . . . , fat) and K = (ki, . . . , k^). A 
pair K, K' will be called allowed if both sets contain N different vectors and 
for any i and j, kj -|- k^ is in A* = (27r/L)Z'^, where denotes the set of 
(i-dimensional integers. We have the following. 

Lemma 1 A pair K, K' is allowed if and only if both sets contain N dif- 
ferent vectors, and there exist d- dimensional integers Ui, and a (complex) 
vector ^ such that kj = (27r/L)nj -|- ^ and k^ = (27r/L)n^ — ^ for every i. 

If both K and K' contain N different vectors of the above form, the pair 
is clearly allowed. On the other hand, if the pair K, K' is allowed, choose 
e.g. ^ = ki. Then k^ + ^ = ki + k^ is in A* for every i and kj — ^ = 
(kj -|- k'^) — + ki) is also in the Abelian group A* for every i. 
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Let us consider the product of two Slater determinants, 
V'k.kKR') = det[e*''^'"]det[e*''^'--] = ^ (_)-(_)'r'e^Ef=i(k.o)+k;,(,pr, _ 

7r,7r'eSjv 

(9) 

For any allowed pair K, K', ipK,K' has the properties (Pl-3). It is also an 
eigenstate of — ^ Aj, provided that kj ■ k^- = for any i,j, a condition which 
can be satisfied above one dimension: Since Xli(^^7r(i) + ^l-'(j))^ ~ + 
k^/^-i(j))^, the [Niy plane waves in (9) may belong to at most A^! different 
eigenvalues. Thus, ^k,k' is an eigenstate if ^^(kj + k^^.))^ = + k^)^ 

or, equivalently, 

N N 

Y.^^■Kii)=J2^^■K (10) 
i=l i=l 

for all permutations tt G Sn, which trivially holds if every kj is orthogonal 
to every k^ . The following theorem suggests that there can be a huge redun- 
dancy in these equations. It extends the above example to further eigenstates 
and shows that for K and K' it suffices to satisfy only (A^ — 1)^ equations. 

Below we use the notation {iii2 . . . im) for the cyclic permutation carrying 
ii into 12-1 etc., im into ii. 

Theorem 1 Let K, K' he an allowed pair. The following assertions are 
equivalent. 

(i) Equation (10) holds for n = (j j+1) with I < j < N—1, forn = {j j+lj+2) 
and (j j+2 with 1 < j < N -2 and for tt = {jlj+l l+l) and {j /+! j+l /) 
with 1< j + 1< I < N. 
(ii) 

(k,+i - k,) • (k;+i - ko = , i<j,i<N-i. (11) 

(Hi) Equation (10) holds true for all permutations. As a consequence, V'k.k' 
is an eigenstate of the kinetic energy belonging to the eigenvalue X^C^i + '^D^ 
whose degeneracy is at least (A^!)^.- '?/'k,k' is the linear combination of {N\y 
different, pairwise orthogonal, N -particle plane waves. Furthermore, each 
occurring plane wave is associated with N different wave vectors chosen from 
{kj + k^}^=i whose all the elements are different. 

Since (iii) obviously implies (i), it suffices to show that (ii) follows from 
(i) and implies (iii). 

(i)^(ii). Observe that we impose (10) for {N — 1)^ different permutations, 
and this is just the number of equations (11) we want to prove. Writing 
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(10) for TT = (j j + we immediately obtain (k-j+i — k^) • (k^-^;^ — k^) = 0. 
If denotes the equation (10) then £'(j^i^2) ~ ^Ur^i) ~ -E'(>+-i>f2) yields 
(ki+2 - kj+i) ■ (k;-+i - k;.) = and Ey^2j+i) - ^o>hi) - Ey+ij+2) yields 
(kj+i - kj) • (k^_^2 - k^+i) = 0. Moreover, E^jij+n+i) - E(^ji) - £^(j4.u+i) is 
equivalent with (k/+i — k;) ■ {^'j+i ~ k^) =0 from which the last equation 
follows by interchanging j and 

(ii)^(iii). First, we note that (11) implies (kj — k|) ■ (k'^^ — k'J = 
for all j,l,m,n. If vr = (j/), (10) is equivalent with (kj — k;) ■ (k^ — kj) = 
0. For a general tt different from the identity we write down the equation 
(10) for TT and for Tr~^ (which may be the same) and rearrange them as 
Eki • {K{i) - ki) = and ■ (k- - k^(i)) = 0, respectively. The sum 

of them yields X^(kj — k^(i)) ■ (k'^^^,.^ — k^ = which holds because each term 
separately vanishes. The difference gives XI (kj + ^Tr^i) ) ■ {^'^{i) ~ k'J = from 
which we subtract 2ki-^(k'^j.^-j — k'J = 0. Again, the resulting equation holds 
because each term equals zero. It remains to show that in ■ip-K,K' there are 
(iV!)^ different plane waves and in each of them N different wave vectors. For 
this we note that all the N'^ wave vectors kj+k^ are different: kj+k^- = k^+k^ 
with 7^ (/, m) would be in conflict with (11). Now the coincidence of two 
plane waves would mean k7ri(j) + k'^'(i) = ^712(1) + ^w^i) ^1 7^ ^2. 'n'l 7^ 
and for all i, which is impossible; nor can we have kj + k^^^^ = k^ + k^j.^.^ 
unless i — j. This ends the proof. 

Let us remark that (11) is equivalent with 

6(k, + k;) + e(k,+i + k;+,) = 6(k, + k;+,) + e{k,^, + k;) (12) 

where e(p) is the one-particle energy, p^. That is, the two-particle energy is 
conserved during an elastic collision of two particles, see also later. 

The general solution of (11) can easily be obtained. The vectors kj+i — kj 
and kj_^i — k'j must be in different, orthogonal, subspaces of M.^. One vector 
from each set, for example ki and k'^, can be arbitrarily chosen from A*. In 
two dimensions both subspaces are one-dimensional and we have 

27r 27r 
kj = ki + —Piik, mo) , k^ = k; + —Piimo, -lo) {i^2,...,N) . (13) 

Here /q and mo are integers at least one of which is different from zero, 
P2, ■ ■ ■ ,Pn are different nonzero integers and the same holds for p'2, . . . , p'^. In 
three dimensions one of the subspaces is one-, the other is two-dimensional. 
Without restricting generality we may suppose that kj — kj are in a one- 
dimensional subspace. Let no be any nonzero integer vector. Since the 
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subspace orthogonal to no has an infinite intersection with Z^, there is no 
problem to choose the vectors of K'. So we have 



ki = ki + — Pino , k; = k; + —n[ (i = 2, . . . , TV) , 



(14) 



where pi are different nonzero integers and n^ are different nonzero integer 
vectors orthogonal to no. 

A complex vector ^ (the same one which has appeared in Lemma 1) can 
be added to ki and subtracted from k'^^ without changing the eigenvalue 
or the eigenfunction: it is a kind of gauge variable which affects only the 
decomposition of V'k.k' in the product of two determinants. 

We note that the two examples given in the Introduction correspond 
to Vk,k' with K = ±^(n,0), K' = ±^(0,m) and K = ±f(n,n), K' = 
±^(n, —n), respectively. 

The question naturally arises whether with Theorem 1 we have described 
all the eigenstates of the kinetic energy operator satisfying (Pl-3). For the 
moment we cannot even exclude the existence of other eigenstates V'k,K' 
which do not fall under the authority of Theorem 1. For a general '?/'k,k' with 
an allowed pair K, K', some elements of {kj + k^}j,j may coincide and, thus, 
the number of different wave vectors can be less than iV^. This may also im- 
ply that, by cancellation, the number of different A'"-particle plane waves con- 
tributing to V'k.k' is less than (Nl^. (Cancellation occurs if {ki + li'^^^^-^yfL-^ — 
{ki + k.'g(^^^}^i for two permutations a and f3 of different parity. It is easily 
verified that for N = 2 and 3 this is not possible but can take place for > 4. 
An example for = 2 and A^ = 4 is K = ^((0, -1), (1, 0), (0, 1), (-1, 0)) and 



K'= ^((-2, !),(!, !),(-!, 0), (0,0)) with a = (24) and /? = (12)(34).) As 



a consequence, some of the equations (10) could not be satisfied but ■0k,k' 
still could be an eigenstate. In the following, we exclude this and all other 
possibilities. 

Theorem 2 Ford > 1. am,ong the eigenstates of the kinetic energy which 
satisfy the conditions (Pl-3) the set of those described in Theorem 1 contains 
a basis. In one dimension such eigenstates do not exist. 

To begin the proof, we note that an eigenstate of the kinetic energy 
satisfying (PI) and (P2) is a finite linear combination of symmetrized plane 
waves. 




(15) 
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In this expression a goes over some index set 7, P'' = (p", . . . , p^), 

^+4R) = ^ e'^^P"«"^- (16) 

7re5;v 

and, for each a, aU p" arc in A* and ^j(p")^ is independent of a. If 
are complex, = V'Re + ^V'lm where ^Re = ZlQe/(^®"«)^P" V'lm = 
^^gj-(Imaa)'?/'pa. It is easily seen that (P3) holds for ip if and only if it 
holds both for -(/^rc and ipim- Therefore, when looking for a basis among 
(Pl-3) eigenstates, we can limit our considerations to functions (15) with 
real coefficients. To further specialize our choice of a basis, we need the 
notion of minimality. When writing (15), we understand that all the possible 
contractions in ip have been made, so that Oq if a G / and the sets 
{p"}^i are different for different a (but may overlap with each other). We 
call the orthogonal set {ip^a}aei of (Pl"2) eigenstates, all belonging to the 
same energy, the support of tj) and denote it by supp^. A (Pl-3) eigenstate 
will be called a minimal function if it is the unique (Pl-3) eigenfunction 
(apart from a constant multiplier) in the linear span of its own support. 

It will suffice to study only minimal functions because of the following. 

Lemma 2 Among the (Pl-3) eigenstates there exists a basis formed by 
minimal functions. 

Suppose that is a non-minimal (Pl-3)-eigenfunction. It is enough to 
show that it can be decomposed as ^ = a0i-|-fe02 where 0i and 02 are (Pl-3)- 
eigenfunctions whose support is strictly smaller than that of 0. Then, if (pi is 
not minimal, it can be decomposed in the same manner as and, in a finite 
number of steps, we arrive at an expression of as a sum of minimal functions. 
Now, because = X^Q,g/ Oa^pa is not minimal, its support spans an at least 
two-dimensional subspace of (P3)-eigenstates. Take here a 0' = X^ae/ ^'a'0p« 
which is linearly independent of (and may have vanishing coefficients). 
Choose a (3 such that 6^ 7^ and define — (j) — ^<f>'- Take a 7 7^ /3 and 

define 02 = 0' — ^0. Then 

0=(l-^4l)""0,+ (^-M"^, (17) 

is the decomposition we were looking for. 

Prom now on, we shall suppose that ip also satisfies (P3) and is minimal. 

We are going to prove that the minimal functions are precisely the eigenstates 
described in Theorem 1. 
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First we consider N — 2. The simplest (Pl-3)-eigenstates are of the form 



^(ri,r2) 



,i(piri+p2r2) _|_ „i{p2ri+pir2) _ „i(qiri+q2r2) 



,i(q2ri+qir2) 



^p(ri,r2) -t/'5(ri,r2) 



(18) 



where Pi, P2, Qi, q2 e A* and 



Pi + P2 = qi + q2 



2,2 2,2 

Pi + P2 = qi + q2 



(19) 



or equivalently, 



Pi + P2 = qi + q2 



Pi • P2 = qi • q2 ■ 



(20) 



These eigenstatcs are certainly minimal, because a (P3)-eigenstate must con- 
tain at least two symmetrized plane waves. 

In any dimension equations (19) have trivial solutions, pi = qi and p2 = 
q2 with pj e A* arbitrary, which are useless because they yield ip{ri, ¥2) = 0. 
Thus, we must impose {pi, P2} H {qi, q2} = 0, since an overlap would imply 
coincidence of the two sets. Moreover, if pi = p2 = p then qi = q2 = 
p, a solution already discarded. Hence, in any nontrivial solution of (19), 
Pi) P2, qi and q2 are four different vectors. 

In one dimension, if any of the four numbers is zero, we obtain a trivial 
solution. In the opposite case q2 = PiP2/<li- This we insert in the first of Eqs. 
(19) to find {pi — qi){qi — P2) = which, again, leads to a trivial solution. 

For d > 1 let pi, P2, qi and q2 be any nontrivial solution of Eqs. (20). 
Choose an arbitrary vector k2 and introduce ki = qi — k.2, k2 = p2 — kg 
and \s.[ = pi — qi + kg. Then ki 7^ k2 and k'^ 7^ kg. The original vectors 
can be written in the form pi = ki + k'^, p2 = k2 + kg, qi = ki + kg and, 
because of the first of Eqs. (20), q2 = k2 + k^ Furthermore, (20) implies 
(k2 — ki) • (kg — k^) = (p2 — qi) • (qi — pi) = 0. Thus, the eigenstate (18) 
takes on the form prescribed in Theorem 1. Clearly, kg plays the role of the 
gauge vector ^. 

It remains to show that minimal eigenfunctions, being the sum of more 
than two symmetrized plane waves, do not exist for N = 2. Let us sup- 
pose the opposite, i.e. that for some n > 2 there is a minimal eigenfunction 



^(ri,r2) = E!^^=i«m^pm, where = {pT,P^) and (p^^)^ + (p™)^ is in- 



dependent of m. Because of the minimality, identical vanishing of ip{r,r) — 
2 Ylm=i flmC^'-Pi''''^™^'" implies that also p™ + p^ is independent of m (indeed 
possible for n > 2, see Lemma 3 below), and Ylm=i'^rn — 0. But then 
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'0(ri,r2) = J2m=i «m['0pm(ri,r2) - '0pn(ri, r2)], where each term of the sum 
is a minimal eigenfunction of the type (18), so ijj could not be minimal. Thus 
we have shown the theorem for N = 2. 

We complete the study of the case = 2 by answering the following 
question: Given pi and P2 7^ Pi, both in A*, how to obtain all the solutions 
of Eqs. (19) for qi and q2? 

Lemma 3 Let pi 7^ p2,. both in A* , and consider the sphere S in 
drawn over pi — p2 as a diameter. There is a one-to-one correspondence 
between nonzero orthogonal pairs of vectors in the intersection of S with A* 
and the solutions of equation (19) for {qi, q2}. //a and b form such a pair, 
the corresponding solution is = pi — a = p2 + b and q2 = P2 + a = pi — b. 

If there is a pair {qi,q2} 7^ {Pi;P2} which solves (20), we use the 
parametrization with an allowed pair K, K' described above and define a = 
k'j^ — kg and b = ki — k2. Then qi and q2 are as asserted, a • b = and 
Pi — P2 = a + b. As we have seen in Lemma 1, both a and b are in A* and, 
because they are orthogonal and sum up to pi — P2, they are also on the 
sphere 

S^{veM.':\v- i(pi - P2)| = ^|pi - P2I} . (21) 

Oppositely, let a e A* fl 5" and define qi = pi — a and q2 = P2 + a. Then 
qi + q2 = Pi + P2 and qi • q2 = Pi • P2 + a • (pi - P2) - a^ = pi ■ p2 because 
a G 5*. Now along with a, b = pi — p2 — a is also in the intersection and 
a ■ b = 0. So vectors in A* n 5' occur in orthogonal pairs whose elements add 
up to pi — P2 and provide the same solution. 

Any nontrivial partition of (27r/L)~^|pi — P2P into the sum of the square 
of 2d integers gives rise to at least one solution of (19). Therefore, if pi — p2 
has two nonvanishing coordinates, there exists at least one orthogonal pair of 
vectors described in the lemma. On the other hand, if only a single coordinate 
is nonvanishing, there may not be any nontrivial solution for qi, q2 (take, e.g. 
d = 2 and pi - p2 = (27r/L)(0, 3)). 

We turn to the case N > 2. Since i/j(R.) is a symmetric function, for (P3) 
to hold ijj{r, r, r^, . . .) = suffices. To see how -0(1", i", i"3, ■ ■ ■) can vanish, we 
write 

^p«(R)= Yl ^p"W (22) 

l<J<m<iV 

with 

■0^(R) = |^g»(P!'"l+Pmr2) _|_ gi(Pmri+Pir2)j g«E^=3P7r(n)''n 

7re5iv:7r(l)=Z,7r(2)=m 
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= ^;p.P.}(ri'r2)^;p.w(r3,...,r^). (23) 

Here we have used again the notation (16) for symmetrized plane waves. 
Now ■^^(r, r, r3, . . .) can be cancelled as a whole by an identically equal 
■0p^(r, r, r3, . . .). The condition of an identical equality imposes 

{Pj}.-^/,m = {P^h^s,t (24) 

and 

pr+p^=pf+pf , (prf+(p^f = (pfr+(pff (25) 

where we have added the last equality which comes from the fact that ip^g 
must belong to the same energy eigenvalue as ippa- So for the vectors 
p",p"^,pf,pf we find precisely the equations (19). Again, we have to dis- 
card the trivial solutions. They would yield P" = P^, so ■?/'pa(R) = ■?/'p^(R), 
which contradicts the supposed contracted form of ijj. Considering equation 
(25) for every pair Z < m, we conclude that for each a the set {p"}^i contains 
N different vectors. Furthermore, if ijjpT'^{r, r, r^, . . .) = ■^p^^ (r, r, r^, . . .) 

then {p'j}f=i n {Pj}f=i = {P'j}j^h,mi = {Pj}j^suti, and for the same a and 
f3 no other choice of Z < m and s < t can yield the equality (24). Thus, 
to annihilate ippair, r, rs, . . .) we need one term from the decomposition (22) 
of N{N - l)/2 different Solving the N{N - l)/2 equations (25) for 

{p'f,pf}, one for each pair / < m, we may (but for the moment we cannot 
say we really do) obtain N{N — 1) different vectors which, together with 
the N vectors of P"", form a set of iV^ wave vectors. Recall that V'k,k' of 
Theorem 1 are also constructed with the help of N'^ vectors kj + k^-. So far 
we have found that minimal eigenfunctions share some attributes of V'KjK', 
but have not seen as yet that these latter themselves are minimal. 

Lemma 4 '0k,k' of Theorem 1 are minimal (PI -3) -eigenfunctions. 

In the present case the index set is Sn- Let (K + a'K')j — kj + k^_i^^.^. 
We have to prove that 

^K,K' = E (-)>K+aK' (26) 

is the unique (Pl-3)-eigenstate that can be obtained as a linear combination 
of the functions '?/'K+aK' with a going over S^. According to equations (15), 
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(22) and (23), a general function over the above basis reads 



{(K+aK'),b=^ 



,+ 



(ri,r2) 




{(K+aK'),b^i, 



(rs, . . .) • (27) 



We may suppose that the coefficient of the identity is nonzero - this can 
always be achieved by renumbering the vectors of K', if necessary - and 
normalize ip so as to have Oid. = 1. In (27) we have regrouped the terms 
which can cancel each other if ri = r2. The point is that, given {a, I < m), 
there is a unique (/3, s < t) such that equations (24) and (25) hold true, 
namely, j3 = {lm)a, s = I and t = m. This is because kj + k^- uniquely 
determines (i, j). Thus, ii ijj is a (P3)-eigenstate then for every a and every 
pair / < m, a(^im)a = —eta- Starting with (/m) G Sn, o,{im) = —1, and 
the products of any two, three, etc. of them appear with coefficients (— )^, 
(— )^, etc., respectively. Since the inversions generate Sn, we conclude that 
Qa = {—)°' for each a e and, thus, ip — ■0k,k', proving that the latter is 
minimal. 

The last step is to show that no other type of minimal eigenf unctions 
exists. This we have seen for two particles and will prove for iV > 3 by in- 
duction. Suppose the claim is verified for any n < N. Let V — ^aei{Pi}iLi 
i.e. the full set of wave vectors used to define the minimal (Pl-3) eigenfunc- 
tion (15). After resummation this function reads 



Here every (pp is a uniquely determined (Pl-3) eigenfunction of — ^j-^ Aj 
which, by the induction hypothesis, is a linear combination oi N — 1-variable 
eigenfunctions of the type described in Theorem 1. (That 0p is actually 
minimal and, thus, equals a single function of this type, will follow from the 
proof.) The sum (28) can be considered as an embedding of any given of 
its terms into a minimal (Pl-3) eigenfunction. The {N — 1)^ wave vectors 
appearing in each minimal component of 0p are of the form kj -|- k^- with 
i < hj < N — 1 and satisfy the orthogonality relations (11). Consider the 
function 




(28) 



N-l 



J2 i-T'"' exp{ipr;v + i J2 (^^(J) + K'{j))r3} 
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in the term p of (28). The symmetrical embedding of the plane wave be- 
longing to (tt, tt') implies the presence, in the expansion of ip, N — 1 other 
plane waves, one for each j, in which p and k7r(j) + k'^,Q-) arc interchanged. So 
k^(j) + k'7r'(j) £ ^ for j = 1, . . . , — 1 and embedding of the {N — l)-particle 
plane wave exp{iprj + E5j^(k^(0 + K'(i))^i} 0k,(,)+k;,^^.j (ri, . . . , r^-i) 
implies that p = k + k', and {(k,k')} U {(kj, k^)}jt"/ satisfy (11). We con- 
clude that in V there is a subset parametrized by an allowed pair K, K', both 
members containing N vectors such that (11) holds, and all the correspond- 
ing symmetrized plane waves ipK+aK' occur in the support of ip. Since ip is 
minimal, this is the support of ip, and -0 = V'k.k'- This ends the proof of the 
theorem. 

2.2 Remarks 

1. The representation of a bosonic state as the product of two Slater de- 
terminants is, in a certain sense, natural. Intuitively, we can imagine a 
spinless boson to be built up of two coinciding half-spin fermions of opposite 
spins, and the two Slater determinants as wavef unctions of the two kinds 
of fermions. Also, the eigenstates we are looking for are twice continuously 
differentiable (C^, in fact, analytic) functions of all coordinates and vanish, 
therefore, at least quadratically at rj = r^. The C^-functions of all coor- 
dinates satisfying (Pl-3) form a linear space C, and one may suspect that 
■0K,K' with allowed pairs is a generating system of this space. A way to prove 
it would be to take an arbitrary function of C, divide and re-multiply it 
with a 'nice' antisymmetric function and expand both - now antisymmetric 
- factors as a linear combination of Slater determinants which are known to 
form a basis in the antisymmetric subspace. The crux is the choice of the 
auxiliary antisymmetric function. It should yield an at least L^-convergent 
expansion for both factors, and the product of the two expansions should still 
be L^-convergent. Because C is dense in the symmetric subspace of L^(A^), 
a successful proof would imply that V'k.k' with allowed pairs linearly span 
the symmetric subspace. 

2. The form of ■0k,k', corresponding to (28), is 

N 

V'K,K'(R) = E (-)'''"'"™e^^'^'^^-^-''"</'k,+k;.(ri, ...,r^-i) , (29) 

l,m=l 

where (f>k^+k>^ = V'K\{kj,K'\{k;,}- 
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3. Theorem 2 is not claiming that eigenstates, being, for instance, the 
product of more than two Slater determinants, do not exist. It is easily seen 
that in d > n dimensions, with n pairwise orthogonal k-sets one can obtain 
eigenstates in the form of a product of n Slater determinants. However, for 
n even these states can be expanded in the basis {^k,k'} (as for n odd they 
can be expanded in the basis of Slater determinants). As an example, the 
expansion of a product of four determinants belonging to K, K^, and K^, 
respectively, reads 

■0K,Ki,K2,K3 = ^ (-)"^'0K,Ki+aK2+^K3 , (30) 

a,l3eSN 

irrespective of the space dimension or whether the product is an eigenfunc- 
tion. 

4. The energy of the eigenstates ipK,^' is easy to estimate. In two and 
three dimensions the length of the vectors in one of the two sets has to 
increase at least linearly, so the energy is at least of the order of N^/L"^ 
(because of the orthogonality (11), at the same time 

^k, • k; = ki • J2iK - K) + K ■ - ki) + Nk, . K 

is only of order N'^/L'^). At fixed positive densities this is much larger than 
the ground state energy of any physically meaningful Hjsf (which is smaller 
than a multiple of N), and these states play no role in thermodynamics either. 
However, in cases when the grand-canonical ground state is the vacuum, 
A*" = 0, as e.g. for repulsive pair interactions and zero or negative chemical 
potentials, we obtain gapless excitations in the form oi N — o(L^/'^)-particle 
eigenstates whose energy vanishes with the increasing volume. 

5. Lemma 3 shows that degeneracy of minimal states can occur. Because 
of nonvanishing overlaps, within a degenerate subspace the minimal states 
may not be linearly independent and one may not choose an orthogonal basis 
of them. 

6. ■0K,K' are eigenstates of the total momentum operator with eigenvalue 
A"q = ^ (kj + k ■) . Multiplying and dividing ■0k,k' by exp i q X] r ^ yields 

^K,K'(R) = e''^^^^^K-k,K'-F(R) (31) 

with k = A^~^ ^ kj and K — k = (kj — k)^^ and similar for the primed 
variables. If V'k.k' is an eigenstate of the kinetic energy then V'K-k k'-F 
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a zero- momentum eigenstate whose k-sets satisfy not only (11) but also the 
more specific equations (kj — k) ■ (k^ — k') =0 which are easy to verify. 
Replacing kj — k by kj and k^ — k' by k^ we obtain the canonical form of 
minimal functions: 

Corollciry A basis of (Pl-3) eigenstates can be chosen among the func- 
tions 

V'q,K,K'(R) = e^'>^'"^V'K,K'(R) (32) 

where q e A*, A^q is the total momentum and ■0k,k' is a zero-momentum 
minimal eigenfunction. Namely, K, K' is an allowed pair with kj • k^. = for 
any i and j and ^ kj = ^ kj = 0. 

In the right-hand side of equation (32) R can be replaced by (rj — r)^;^ 

where r = Z^ili ''j- 

7. Equations (19) and (25) suggest to interprete the (Pl-3) eigenstates of 
the kinetic energy operator as stationary states of a system of impenetrable, 
pointlike bosons which interact with each other through elastic collisions. 
Clearly, there exists no self-adjoint Hamiltonian with a dense domain in the 
symmetric subspace of L^(A^) which would describe such a system! The 
states we have found cover only a very small fraction (roughly ~ (A"!)~^) 
of the Hilbert space. Since the suitable operator would be Hi\f with infi- 
nite repulsive delta-interactions, Theorem 2 provides an independent proof 
of the otherwise known fact, cf. the Introduction, that above one dimension 
— ^Aj, defined as a symmetric operator on C^(A^ \ Uj^j{rj = r^}), 
has no self-adjoint extension describing such an interaction. The states char- 
acterized in the above theorems are in the domain of H^, the closure of 
H^. In two and three dimensions they are eigenstates of any Hamiltonian 
H]sf with a properly defined contact interaction, and the only eigenstates in 
the domain of Any such Hn is a self-adjoint extension of or H'^, 
so densely defined in L^(A^), and having infinitely many other symmetric 
eigenstates, which do not vanish (even diverge) at particle encounters. As 
suggested in the Introduction, we can orthogonally decompose any Hf^ in an 
interaction-dependent part and i^eiastio cf. equation (7). 

8. Because of their simple form, some basic properties of minimal func- 
tions can be studied easily. Nodal surfaces, for example, can be identified if 
the function is written in the canonical form (32). Let V and V be the or- 
thogonal subspaces containing kj and k^-, respectively. Since dct[e*'^''^'"] = 
if Fj — Tj is in V' and detfe*''''''"] = if rj — Vj is in V, both for any i ^ j, 
we find that t/'K,K'(R) = if r j — is in V U V for a pair i 7^ j. So there is 
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a regular long-range exclusion effect, although it concerns only a set of zero 
Lebesgue measure in the configuration space. This effect can be amplified 
by particular choices of K or K', that we can see easier on the two-point 
correlation function. A straightforward computation of this latter yields 



P2(0,r) = 7V(Ar-l)(5(ri)5(r2-r)) 



N 



N-1 



1 




,ik{r 



(33) 



where p — N/\A\ is the uniform value of pi(r) = N{5{ri — r)) and av- 
eraging means integration over ri,...,rAr in A with the weight function 
|V'k,k'|^/||'V'k,k'||^- III the derivation of (33) the minimality of V'k.k' is ex- 



m. 



ploited via the property that kj -|- = k; -|- k^ implies i — I and j 
For most choices of K and K', p2(0,r) outside a small neighborhood 

oi V U V where it vanishes. In general, let be the set of points r in A 
such that k; ■ r is an integer multiple of 27r for every /. is an Abelian 
group with respect to addition modulo A which contains V as a subgroup 
and may contain translates of V. For example, for i = 1, . . . , o? let Mj be the 
greatest common divisor of {(L/27r)/c„j}^^jwhere kni is the ith component of 
k„. Then contains the points with coordinates Xj = rriiL/Mi, where rrii 
is any integer between —Mi/2 and Mi/2. Defining Qj^i analogously, we get 
another Abelian group containing V and maybe some translates of it. Now 
P2(0, r) = if and only if r is in U ^k' and, as we can check also directly, 
V'K,K'(ri) = if (and only if) r, — is in U ^k' for a pair i ^ j. 

9. In the cigcnstatcs ^/'k.k' there is no off-diagonal long-range order: 
The associated one-particle reduced density matrix, a, has A^^ nonvanishing 
eigenvalues, each of which is equal to 1/A^. Indeed, using that all k; + k^ 
are different, a straightforward calculation of the integral kernel of a yields 



<7(r,r') 



^iiV'K,K'ir 

N A ^ 

Lm=l 



' / V'K,K'(r,r2, . . .)V'K,K'(r',r2, ■ ■ ■)dr2 . . .dr^r 



i(k,-hk;,)(r-r') 



(34) 



Because ■0k,k' is a momentum eigenstate as well, a is diagonal in momen- 
tum representation and its eigenvalues are the diagonal elements. Fourier 
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transforming a{r,r'), 



c^(P, P) = ^ ^PMi+i^'m ■ (35) 

l,m=l 

Since p = k; + can hold for at most a single pair {l,m), the result follows. 

10. The theorems extend to systems in Li x • • • x rectangles with 
periodic boundary conditions. The appropriate definition of A* is (27r/Li)Zx 
• ■ ■ X (27r/Ld)Z. The orthogonality relation (11) holds for example if {kj — 
kj} and [k'^ — k'j} have disjoint subsets of nonvanishing coordinates. The 
existence of other solutions depends on the rationality of the ratios among 
Li. 



3 Eigenstates for lattice models 

We consider the same problem as before on an L x ■ ■ • x L part (now L is 
a positive integer) of Z'^ with periodic boundary conditions. Defining the 
lattice Laplacian as 

{A^){v)^ J2 mr')-m], (36) 

|r'-r|=l 

e^^^ with k — (ki, ... ,kd) is an eigenfunction of —A belonging to the eigen- 
value e(k) = 2 X]f=i(l "COS ki). A pair K, K' is called allowed if both {k^}^^ 
and {k'j}jLi contain N incongruent vectors modulo 27r and kj + k^- is in A* for 
all Now for an allowed pair K, K', ipj^ k' as given by (9) is an eigenstate 

of - Ef=i A. if 

N N 

E^(ki + k;«) = E^(ki + kD (37) 

i=l i=l 

for all TT G Sn, in which case (37) provides the eigenvalue. The analog 
of Theorem 1 can be derived, based on the following observation. For 2d- 
dimensional complex vectors u and v let us introduce the dot product u • v = 
Y^^i ^^ciA, is, without taking the complex conjugate in one of the factors. 
Then, setting 

v(k) = (cos /ci, . . . , cos /crf, i sin /ci, . . . , i sin kd) , (38) 
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the condition (37) is equivalent with 

N N 

5]v(kO-v(k;(,)) = 5]v(k,).v(kD. (39) 

i=l i=l 

The analog of the particular solution kj • k^ = in the continuum case has 

to satify v(kj) ■ v(k^) = or YM=i^^^i^ii + ^'ji) — hJ- For L odd 

no allowed pair can solve these equations. For L even, in two dimensions, a 
family of solutions is obtained by choosing ki2 = 7r/2 — kii and A;', = 71/2 — k^^. 
If 4 is a divisor of L, the same choice works in three dimensions if ki^ + /c^-g = 
7r/2 for all e.g. fcja — fc-g = 7r/4 for all i. More solutions, in particular, 
solutions also for L odd are provided by the following theorem. 

Theorem 3 Let K' he an allowed pair. The following assertions are 
equivalent. 

(i) Equation (39) holds for the permutations listed in part (i) of Theorem 1. 
(n) 

[v(k,+i) - v(k,)] • [v(k;+,) - v(k;)] = , 1 < z, J < AT - 1 . (40) 

(Hi) Equation (39) holds for all permutations. As a consequence, ■0k,k' is 
a (Pl-3) eigenfunction of the kinetic energy and, thus, of all Bose-Hubbard 
Hamiltonians with purely on-site interactions. The eigenvalue is given by 
equation (37). 

The proof is exactly the same as that of Theorem 1; it suffices to replace 
everywhere kj by v(kj) and k^ by v(k9. Also, condition (40) expresses energy 
conservation in two-particle collisions, and is the same as equation (12), if 
we use the suitable expression for the one-particle energy. 

To find solutions of (40), it is convenient to rewrite it as 

d 

^[cos{ki+ii + k'j+ii) + cos{kii + k'ji) - cos{ki+ii + k'ji) - cos{kii + /c^+ij] = . 
1=1 

(41) 

We will consider two families of solutions. Type 1 solutions exist for any L. 
In two dimensions they are obtained by choosing an allowed pair with ka 
and A;^2 independent of 1 In three dimensions we can choose e.g. ka, ki2 and 
independent of i. In these examples equation (41) holds by a separate 
cancellation of the d terms. Also, notice that the pair K, K' satisfies the 
equations (11). Rewriting the solution in the canonical form (31) or (32), the 
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constant components of the transformed vectors vanish and kj • k^- = for all 
Recalling remark 8 of the former section, this implies that '?/'k,k'(-R') — 
whenever for some i,j the vector Tj — Tj is parallel to a coordinate axis. Thus 
we obtain: 

Corollciry Type 1 solutions of equation (40) are eigenstates of the more 
general lattice Hamiltonian (8). 

Type 2 solutions exist only for L even and are given by ki2 = d — kn and 
K2 = TT — -i? — which, in three dimesions, is completed with e.g. ki^ chosen 
to be independent of i. Since the choice of 'd does not influence the solution, 
we can fix 1? = 7r/2, as earlier. 

The solutions of (41) we have presented above exist only for N < L. In- 
deed, in all these examples either in K or in K' a single component specifies 
the vector. Since in both sets all the vectors must be incongruent, the defin- 
ing component must take on different values out of L possible ones. In 
higher dimensions the same kind of solutions can easily be given - for exam- 
ple, in four dimensions we can fix the first two components in the K set and 
the second two in the K' set. In general, we can freely choose N incongruent 
vectors from a [d/2] dimensional subspace, which limits the number of par- 
ticles to < LW'^y If iV -C LW'^\ it is possible to choose K, K' such that 
|ki|, |k^| -C 1 for each i. Then, in (37) we can expand the cosine functions 
up to second order, and find that the smallest attainable energy is of order 

, . 2 2[d/2] 

L~^N '^W2], Thus, for N — o(L'^+W'^]) there are eigenstates whose energy 
tends to zero as L goes to infinity. They correspond to gapless excitations 
above the grand canonical ground state when this is the vacuum (A'" = or 
parallel spins). 

The reader may notice that part (iii) of Theorem 3 claims less than that 
of Theorem 1. It is because equation (40) does not imply that kj -|- k^- are 
different for different pairs This nevertheless follows from (11) for type 

1 solutions, but a counterexample of type 2 is easily obtained. Let L be even 
and choose ki2 = 7t/2 — kn, {ki^ independent of i if d = 3) and k^ = kj. 
As a special choice of type 2, the pair K, K' = K solves equation (41) and, 
for any (i, j), kj -|- k^- = kj -|- k^. The same example shows that in the sum 
over the A^! symmetrized plane waves yielding the (Pl-3) eigenstate ■0k,k 
(cf. equation (26)) there are coinciding terms: '?/'K_|_aK' ~ "^^-ik+k' i^ipli^s 
that '0K+aK ~ V^K+zSK P — ot~^ ■ Now a differs from its inverse if and only 
if it contains at least one cycle longer than two. Thus, our example works 
for all such a (but we need N > 2) and (5 — a~^. As a permutation and its 
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inverse have the same parity, the two occurrences of the same symmetrized 
plane wave arrive with the same sign in '?/'k,k- 

Keeping the earher definition of minimahty, Lemma 2 remains vahd: Min- 
imal functions form a basis among the (Pl-3) eigenstates of — ^ A,. Type 
1 eigenf unctions are minimal, as it is seen by applying Lemma 4. We believe 
that all solutions of equations (40), and only them, are minimal functions. 
However, we cannot prove this and have only a weaker analog of Theorem 2. 

Theorem 4 If ip is a minimal eigenf unction o/ — ^ Aj then supp ip C 
{V^K+aK'laeSjv /^^^ somc allowed pair 'K^K! satisfying equations (40). 

For N — 2 the assertion is equivalent with the stronger claim of Theorem 
2. The second of equations (19) is to be replaced by 

d d 

^(cospij + cosp2j) = ^(cos qij + cos q2j) ■ (42) 

Again, the trivial solution {pi,P2} = {qi,q2} is to be excluded, but now 
the four vectors need not be different. Indeed, with our earlier example 
ki2 — 7r/2 — kii and k'^ = kj for z = 1, 2, Pi = 2kj and qi = q2 = ki + k2 is 
a nontrivial solution. When parametrizing a solution (18) with an allowed 
pair K, K', equation (42) coincides with (40) written for i = j' = 1. The 
minimality of such a solution is trivial, and the proof that no other minimal 
function exists is the same as in Theorem 2. The weaker assertion for a 
general N can be shown as in Theorem 2 by starting with a weaker induction 
hypothesis. 
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